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Abstract. We consider constraints on inflation driven by a single, minimally coupled 
scalar field in the light of the WMAP5 dataset, as well as ACBAR and the SuperNova 
Legacy Survey. We use the Slow Roll Reconstruction algorithm to derive optimal 
constraints on the inflationary parameter space. The scale dependence in the slope of 
the scalar spectrum permitted by WMAP5 is large enough to lead to viable models 
where the small scale perturbations have a substantial amplitude when extrapolated 
to the end of inflation. We find that excluding parameter values which would cause 
the overproduction of primordial black holes or even the onset of eternal inflation leads 
to potentially significant constraints on the slow roll parameters. Finally, we present a 
more sophisticated approach to including priors based on the total duration of inflation, 
and discuss the resulting restrictions on the inflationary parameter space. 
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1. Introduction 

Cosmology has long had the goal of reconstructing the inflationary potential from 
observational data [JJ [2j [3j H]. This was originally envisaged as a two-step process. 
One would first obtain constraints on spectral parameters such as the scalar and 
tensor spectral indices, n s and n t , their runnings, and their overall amplitudes. Armed 
with this information, one could invert the inflationary formulae for these variables 
to obtain the slow roll parameters. These parameters involve derivatives of the 
potential and would provide a Taylor expansion for the segment of the inflaton potential 
traversed as astrophysical perturbations leave the horizon. It is now apparent that the 
optimal approach to reconstruction is to insert the slow roll variables directly into the 
cosmological parameter set; a process known as Slow Roll Reconstruction [6J [TJ, [8]- 
One thus never computes n s et al, but instead fits directly to the slow roll parameters 
themselves. The flow hierarchy [10l [HJ [12], 13] automatically accounts for the scale 
dependence of the slow roll parameters. Reconstruction thus explicitly realizes the hope 
that the early universe can be used as a "laboratory" for high energy physics: the 
slow roll parameters can be determined by GUT or stringy physics, and we are directly 
constraining their values with our Markov Chains. 

The flow hierarchy can be solved exactly to yield the underlying inflationary 
potential as a function of the first M non-zero slow roll parameters |14j . We can thus 
compute the amount of inflation remaining after a specific mode leaves the horizon, as 
a function of the slow roll parameters. This provides a self-consistency check: fitting 
to M slow roll parameters carries the tacit assumption that these parameters fully 
characterize the portion of the inflationary epoch during which observable modes exit 
the Hubble horizon. Under the assumption that the same parameters characterize the 
inflaton dynamics till the end of inflation, one obtains even stronger constraints by 
requiring a sufficiently long inflationary epoch in order to solve the usual cosmological 
problems. In fact, given the quality of present-day data, much of the power of Slow Roll 
Reconstruction derives from priors based on the number of e-folds (Ej, especially when we 
go beyond the first two slow roll parameters. Further, thanks to the analytical solution of 
the slow roll hierarchy, our one use of the slow roll approximation is the computation of 
the perturbation spectrum. This computation can be performed precisely by numerically 
solving the mode equations, as pointed out in [51 [15] and implemented in [TBI fTT] . 
However, even with the WMAP5 dataset, the quality of current cosmological data is 
not high enough to require the use of an exactly computed spectrum, and we do not 
employ it here, instead using the Stewart-Lyth formulae [18] for the power spectra along 
with the exact solution to the background given by the flow hierarchy. 

There are a number of antecedents to Slow Roll Reconstruction. As noted above, the 
"inverse problem" posed by the extraction the inflationary potential from the observed 
power spectrum was first discussed over 15 years ago. Leach and collaborators [T9l[20l[2T] 
write the spectral indices as functions of the slow roll parameters at a given pivot, 
and then constrain the slow roll parameters with data. If one includes higher order 
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slow roll terms when computing the spectral indices, this approach approximates a full 
solution of the slow roll hierarchy. However, since the full potential is never computed, 
one cannot eliminate sets of slow roll parameters for which the number of e-folds is 
unacceptably small. Conversely, Monte Carlo Reconstruction [TOj CHI [12] specifies an 
inflationary trajectory, and then uses the flow equations to move to a point some 50 or 60 
e-foldings after the observable scale leaves the horizon. However, this process does not 
provide an unambiguous spectrum when inflation apparently continues indefinitely, and 
implicitly rules out the possibility of a hybrid transition in models where the potential 
does naturally lead to the termination of inflation. This ambiguity is discussed in passing 
in [11] . and we will return to it below in Section [H Moreover, as originally implemented, 
Monte Carlo reconstruction did not work directly with likelihood functions derived from 
astrophysical data, and thus did not weight different candidate potentials according to 
their agreement to astrophysical data. 

In this paper, we improve the implementation of Slow Roll Reconstruction by 
developing a physically realistic e-fold prior [8J. This constraint is a strong function 
of the e parameter, which specifies the relative amplitudes of the primordial tensor 
and scalar spectra. Secondly we present constraints on the inflationary parameter space 
derived from the recent CMB (Cosmic Microwave Background) data: the 5- year WMAP 
dataset (WMAP5) [22j [231 [1Q [25] and the 2008 ACBAR release (ACBAR) [26], as well 
as high redshift type la supernovae from the SuperNova Legacy Survey (SNLS) [27J. 
In addition, we show that the constraints on the scale dependence ( "running" ) of the 
scalar spectral index derived from WMAP5 data allow models in which the amplitude of 
the density perturbation approaches unity before the end of inflation. These scenarios 
can be ruled out on physical grounds - either through the overproduction of primordial 
black holes [281 [221 EH EH E2J EH Ell E5], or because the inflaton is moving so slowly 
that its motion would be dominated by its quantum fluctuations, rather than its semi- 
classical rolling. The latter situation is equivalent to eternal inflation [36], and inflation 
cannot terminate coherently over a volume large enough to contain our visible universe. 
The parameter space is tightly constrained by excluding scenarios that lead to the 
overproduction of primordial black holes. Likewise, a substantial region of the parameter 
space which permits the onset of eternal inflation is permitted by the WMAP5 data, 
but most of this region is already excluded by the combination of WMAP5 and SNLS. 
This constraint is more important to the WMAP5 dataset than it would have been with 
WMAP3, which excluded a positive running at the 2a level. The central value for the 
running obtained from analyses of the WMAP5 dataset is still negative, but the overall 
preference for a negative running is significantly diminished. 



The Inflationary Parameter Space after WMAP5 



4 



2. Slow Roll Reconstruction 



2.1. Formalism 

For inflation driven by single, minimally coupled scalar field in a spatially flat FRW 
universe the equations motion can be written with the field as the independent variable 
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where mpi = 1.22 x 10 19 GeV is the Planck mass. The HSR [Hubble slow roll] parameters 
e \n are defined by the infinite hierarchy of differential equations [10] 
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The usual slow roll parameters are i] = 1 Xh and £ 
so that £ Xh = for all £ > M at some 
truncated at order M, the hierarchy can be solved explicitly [H], allowing us to obtain 
an exact expression for the potential. Consequently, we can then determine the value 
of at which inflation ends and the number of e-folds that remain when = 0*. 

This formulation includes the implicit assumption that the scalar field dynamics 
is well-approximated by the truncated hierarchy. We do not insist that the truncated 
hierarchy approaches the exact limit as M — > oo, but only ask that the truncated H{4>) 
is an asymptotic expansion of the exact H(<f)) for values of which are astrophysically 
relevant, with M less than the value at which the asymptotic expansion would begin to 
diverge. This covers a very large class of potentials, but excludes "step" models where 
each successive slow roll parameter is larger than its predecessor in a very small range 
of field values [H]. 

The scalar and tensor perturbation spectra are given by [18J 
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where C = —2 + In 2 + 7 w —0.729 and 7 is the Euler-Mascheroni constant. The scale 
dependence of the spectra then follows from the scale dependence of e, rj and H. In our 
Monte-Carlo Markov Chain (MCMC) analysis, the power spectrum is fixed by setting 
A s = P-ji at a "fiducial" scale k+, where we also specify the initial values of the slow roll 
parameters. We then compute via equation ([2]), which then fixes the overall scale of 
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inflation. In practice A s is effectively fixed by the data, putting tight bounds on H^/e+, 
since e*,?y* <C 1. If e* is very small, if* is also small, lowering P h relative to P-^., and 
reflecting the well-known inflationary consistency condition. We take k* = 0.02 Mpc -1 , 
which was the optimal value for WMAP data found in [6j HJJ . Moreover, a useful feature 
of Slow Roll Reconstruction is that we can post-process our Markov chains using the 
flow hierarchy to move our constraints to an arbitrary k* [6]. 

2.2. The e Parameter and the Duration of Inflation 



Consider the "connection equation" for scales in a universe which inflated, reheated, 
and passed through matter-radiation equality, 

k a^H^ *r/7 \l ^end ^reh Hfc a e nHea /rrx 

exp[-JV(fc)] — — 77 — -^A (7) 



ao-Po a 0-^0 a reh a eq H eq a if() 

Here Hk is the value of the Hubble parameter when the mode with comoving 
wavenumber k leaves the horizon, a subscript refers to the present day value, while "eq" , 
"reh" and "end" denote values at matter-radiation equality, reheating (i.e. the point at 
which the universe thermalizes after inflation) and the end of inflation. Now assume that 
the universe is effectively matter dominated between the end of inflation and reheating, 
so p oc a -3 . Between reheating and matter-radiation equality, the universe is radiation 
dominated and p oc a~ A . Equation ([7]) becomes 



a oHo J 3 VPend/ 4 \Preh 



Using the first Friedmann equation in a flat universe, H oc y/p, along with the usual 
convention do = 1, we rewrite the last term as an expression involving p e q and p C rit- 
Recall that a eq = 4.15 x 10~ 5 /Q m h 2 [35], where fl m is the present matter density, and 
h is the Hubble parameter in units of 100 kms" 1 Mpc -1 , and p e q/p C rit = 2f2 m /ag q , so 

N(k)= -Inf^rV W^l+ ik ^ ^ 
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+ In \^-J + ln 220 n mh. (9) 

Up to rounding issues, this formula is identical to equation (6) of [4_6]. During inflation 
p ~ V((p), but there is little to be gained from this substitution here, and we express 
our free parameters in terms of H. We have 

H cq = 5.29 x 10 6 h 3 Vl 2 m H = 9.24 x 10~ 55 h 4 n 2 m m P1 , (10) 

recalling that H = 1.75 x 10~ 61 h mpi in natural units. Substituting for H eq and p eq , 
the Q m dependence cancels. After reheating, we can specify the energy density of the 
universe in terms of the temperature T, but the relationship between T and H depends 
on the number of degrees of freedom. Numerically, H ~ p 1 / 2 /mpi, so while p 1 / 4 and H 
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have the same units, H is numerically smaller than p 1 / 4 or T. We specify k in Mpc 1 , 
and with c = l,H Q = 100/i/(2.99792 x 10 5 ) Mpc" 1 , so 

N(k) = - In + 5 In (^) - \ In + In f ^ 

+ 59.59. (11) 

In practice, Hk and i? e nd are usually very similar, but reheating need not happen 
efficiently, so if re h can be far smaller than H cn a- If the universe is in thermal equilibrium 
above the neutrino freeze-out temperature (around 1 MeV), the present number density 
of neutrinos is a function of the number density of photons, which is given by the 
absolute temperature of the CMB. By default, the Boltzmann code CAMB j?7] we use 
to compute the power spectrum assumes three (almost) massless neutrino species. Thus, 
to self-consistently compute CMB spectra within our Markov chains, we need T re h ^ 1 
MeV. The neutrino sector makes a subdominant contribution to the temperature 
anisotropies, but nucleosynthesis can easily be disrupted by departures from thermal 
equilibrium when T ~ 0.1 MeV [H]. Given the prefactor of 1/6 on the i? re h term, the 
effective difference between these two scales is very small. Thus while one might take 
T re h ~ 0.1 MeV as a physical lower bound on reheating, it is only slightly more restrictive 
to impose T re h ~ 10 MeV, which also guarantees a thermalized neutrino sector. 

We use this information to eliminate combinations of slow roll parameters that 
cannot produce a sufficient amount of inflation, assuming that the same number of 
parameters used to fit the data provide a description of the inflaton dynamics till the 
end of inflation. As the above analysis shows, this quantity depends on H^, H CIl a, and 
our assumptions about reheating. We could also reject models that produce too much 
inflation - but this assumes that inflation ends via the gradual "erosion" of slow roll, 
rather than via a hybrid transition - a valley or a cliff, in the terminology of [8] 
and we do not make this cut, instead treating equation (TlTj) as a lower bound on N. 
Physically, we are assuming that the equation of state between the end of inflation and 
reheating is somewhere between and 1/3. We are thus ruling out a secondary period 
of accelerated expansion such as thermal inflation [19] or post-inflationary cosmic string 
networks [50] - both of which would reduce N. Conversely, a "kination" dominated 
phase where w ~ 1 would increase N beyond the lower bound we compute |51j . In a 
given model of inflation, N is maximized when if re h = H en d and the universe thermalizes 
instantaneously after inflation. Conversely, if thermalization happens slowly, we still 
need T reh ~ 10 MeV. Here we compute H by assuming p rch ~ Tf eh . This is enough for 
our purposes, since the numerical factor which depends on the number of degrees of 
freedom appears logarithmically. In Figure [1] we show the exclusion regions found for 
the {e, r]} parameter space with three versions of the prior - we will see later that this 
cut is visible in the two parameter {e, rj} chains. 

In what follows, we explore the consequences of four different assumptions about 
the duration of inflation. The first is the bare prior that the inflation lasts long enough 
for modes which contribute to the CMB to leave the horizon - without this we cannot 
self-consistently use the inflationary expressions for the power spectrum, and it amounts 
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Figure 1. The regions of the {e, 77} plane excluded by requiring N(k+) > 15 (black), 
T rc h > 10 TeV (dark) and instantaneous thermalization (light). The other prior used 
in this work, T rc h > 10 MeV, looks very similar to the T re h > 10 TeV prior. For this 
hgure only, the amplitude for the scalar spectrum at the fiducial scale has been set to 
the WMAP5 best fit. 



to demanding that N The next weakest constraint is to demand T re h > 10 MeV, 

which protects nucleosynthesis and guarantees a thermalized neutrino population. A 
stronger assumption is to require T reh > 10 TeV, which ensures reheating occurs well 
above the electroweak scale. Finally, we consider "instant" reheating where i? re h = -f^end- 
Since H decreases strictly with time, the following sequence of inequalities must hold: 

H k > H cnd > H Teh , (12) 

where the first relationship is guaranteed by the dynamics of inflation. If we assign a 
numerical value to T reh , we must enforce the second inequality explicitly. Since we put 
a flat prior on e, our chains typically sample very few points with e < 10~ 5 . The initial 
value of e and the observed amplitude of the spectrum fix H^, so models with low e 
also have a low inflationary scale. For the explicit values of T re h we use, this constraint 
is always satisfied in practice. However, if one puts a logarithmic prior on e, Hk can 
be minute - and any numerical bound on T re h implicitly provides a lower bound on 



I This number is somewhat arbitrary, but it cannot be made significantly lower while simultaneously 
ensuring that the perturbation spectrum is computed self-consistently. 
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Primordial Parameters 


Dataset Combination 


e-fold Prior 


{e^logflO 10 ^} 
{e,?^,log[i0 10 A s ]} 


WMAP5 
WMAP5 + ACBAR 
WMAP5 + SNLS 


A* > 15 
T reh > 10 MeV 
T reh > 10 TeV 

H e nd -^reh 



Table 1. MCMC parameter fits described in this paper. We run chains for all 
combinations of variables, dataset and e-fold prior. 

log(e). When e is very small, it effectively decouples from the slow roll hierarchy [8], so 
parameter regions found with a flat prior on log e would be determined by the prior and 
not the data. Consequently, we do not pursue this question here. 

3. Parameter Sets and Monte Carlo Markov Chains 

We carry out MCMC studies of two sets of "primordial" parameters: {e, rj, log[10 10 A s ]} 
and {e, rj, £, log[10 10 v4 s ]} with flat priors on each of these variables, and the 
standard "late-time" cosmological parameters are also allowed to vary in the 
chains. In the language of [8] these are High-e 2-Parameter and High-e 3- 
Parameter models, respectively. The first parameterization is roughly analogous 
to the standard {r, n s , log[10 10 ^4 s ]} formulation; the second corresponds to the 
extended {r, n s , dn s /d In k, log[10 10 A s ]} set. We constrain these models with three data 
combinations: WMAP 5 year data (WMAP5), WMAP5 combined with the ACBAR 
2008 release (WMAP5 + ACBAR) and WMAP5 combined with the SuperNova Legacy 
Survey (WMAP5 + SNLS) 1§] We do not marginalize over the amplitude of the Sunyaev- 
Zel'dovich (SZ) effect with a flat prior and an SZ template for a fixed set of cosmological 
parameters as was done with WMAP5 in [21]. In principle each vector of cosmological 
parameters visited by the chains has a unique SZ spectrum which should to be added to 
the primordial Ct before being compared with the data, but computing the SZ templates 
is not a trivial task. In practice, for all relevant parameter fits in [23], the posterior of the 
SZ amplitude parameter essentially reproduces the prior, suggesting the SZ effect is not 
detected in WMAP5. Conversely, marginalization over the SZ parameter can potentially 
bias other parameters via volume effects, which depend on the chosen width of the 
prior. For these reasons we choose not to include an SZ parameter in our chains. For 
consistency, we only consider ACBAR bandpowers at £ < 1800, since higher multipoles 
are potentially SZ-contaminated. A summary of our parameter fits is given in Table [U 
We use modified versions of the Boltzmann code CAMB |47] and the public 
COSMOMC code |52j for our MCMC parameter estimation. Each run uses eight 
chains. To assess the convergence of the chains, we apply two stringent requirements: a 
Gelman & Rubin [53] criterion on the least-converged eigenmode of the full parameter 

§ We checked the combination WMAP5+ACBAR+SNLS but found that the ACBAR data did not 
significantly change results relative to WMAP5+SNLS and thus do not present those constraints. 
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covariance, R — l< 0.01, and a conservative criterion on the convergence of confidence 
limits, stopping when the worst rms difference in the 95% CL between chains for all 
parameters is 10% of the standard deviation. 
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Figure 2. The joint 68% (dark) and 95% (light) confidence levels obtained on the first 
two HSR parameters e and rj at the fiducial scale k+ = 0.02 Mpc -1 from the WMAP 
5 year data, assuming £ and all higher order slow roll parameters are zero. The grey 
(dotted) contours correspond to a "minimal" e-fold prior, N(k+) > 15. The blue (solid) 
contours are obtained by assuming instantaneous reheating. 



Our constraints on the inflationary parameters for different data combinations and 
various reheating priors are presented in Tables [ATI and lA2l The WMAP5+SNLS data 
combination is significantly more constraining than WMAP5+ACBAR. We refer the 
reader to these tables for quantitative details, and now highlight the key points. The 
direct effect of the reheating prior on the {e, rj, log[10 10 A s ]} parameterization is a "bite" 
into the constraints, as displayed in Figure [TJ We illustrate the impact of the choice of 
reheating prior on the parameter constraints derived for WMAP5 in Figure [^jj]. In Figure 
El we show the constraints in the {e, rj} plane for a representative (10 TeV) reheating 
prior, for both the WMAP5 and WMAP5+SNLS data combinations. We transform 

|| The keen-eyed reader may notice that the cut in Figure [T] is slightly higher than what appears in the 
constraints in Figure [2] This is an unavoidable artifact of the binning used to calculate 2D contours 
from MCMC using the public GetDist statistics package. The MCMC itself reproduces the sharp cut 
of the prior, but the binning aliases this curve to slightly lower values. 
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these constraints into the empirical {r, n s } parameters by using the approximate second 
order slow roll formulae, and compare these to constraints obtained directly from an 
MCMC analysis where {r,n s } were varied in the chains. The difference follows from 
the choice of e-fold prior: if the prior eliminates regions which would otherwise have 
a high likelihood, the resulting confidence levels will stretch further into the tails of 
the likelihood distribution, where the likelihood is somewhat lower. This can lead to 
a broadening of the marginalized constraints if the slope of the likelihood surface is 
relatively shallow in the tails, and this effect can be seen clearly in the 68% contour 
in Figure [2j We also show trajectories as a function of the number of e-folds for three 
common slow roll models. The A0 4 "toy-model" is firmly excluded, while m 2 </> 2 and 
natural inflation models are still allowed at the 95% CL level for N > 50 - which 
accords with the number of e-folds we would expect for these models. We find that 
our WMAP5+SNLS analysis disfavors a blue spectral index even when r (or e, in our 
case) is included in the chains, which accords with the WMAP Team conclusions for the 
data combination of WMAP5, a compilation of supernovae data and baryon acoustic 
oscillations data [25] . 

Figure H] shows the difference between a minimal e-fold prior and the strongest 
(instantaneous) reheating prior for the {e, r], £, log[10 10 A s ]} parameterization. The e- 
fold priors tend to reject models with positive £, or negative running. Even a minimal 
(N* > 15) prior puts tight restrictions on the allowed upper limit of £. Adding a low 
(> 10 MeV) reheating prior tightens this effect substantially, but adding a stronger 
reheating prior does not affect this limit further, implying that inflation ends so quickly 
for these models that the weakest reheating prior is sufficient to cut them out. We 
will return to this topic in Section HI Figure [5] shows constraints on {e,n,^} for a 
representative (10 TeV) reheating prior for WMAP5 and WMAP5+SNLS. Figure [6] 
shows constraints for the same chains, after being post-processed (maintaining flat priors 
on the HSR parameters) to express them in terms of empirical power law parameters, 
using second-order slow roll formulae. We compare these results to a standard analysis 
using the power law parameters {n s , r, dn s /d In k} in the chains. We see that the e- 
fold prior excludes most of parameter space allowed by naively fitting dn s /d\nk to the 
data. As is now well-known, a large negative running (corresponding to large positive £) 
does not provide sufficient inflation. Marginalizing over this degeneracy with dn s /dlnk 
shifts the constraints on n s and r from the empirical analysis to redder spectral indices 
compared to Slow Roll Reconstruction with a reheating prior. The difference between 
the empirical constraints shown in Figure [6] and those from the WMAP Team's analysis 
in [24l [25] is due to the fact that they used a pivot scale of 0.002 Mpc -1 (compared with 
our fiducial scale 0.02 Mpc -1 ), which exacerbates this degeneracy and hence shifts the 
marginalized constraints through a volume effect. 

The exclusion of large positive £ along this degeneracy by the e-fold prior again 
causes our MCMC to explore parts of likelihood space which were not within the 
marginalized constraints before this prior was applied, and allows a significant negative 
£. We will explore the consequences of the lower bound on £ in the next section. 
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Figure 3. The joint 68% (inner) and 95% (outer) bounds on the slow roll variables 
(top) and power law spectral parameters (bottom), with k+ — 0.02 Mpc -1 , for a High-e 
2-parameter fit. WMAP5 constraints are blue, and WMAP5+SNLS constraints are 
red. Solid contours come from Slow Roll Reconstruction with prior T re h > 10 TeV. The 
running of the scalar index in these models is second order in slow roll and hence very 
small. The dashed contours show results from WMAP5 chains where the spectrum 
was specified via n s and r. We superimpose the "trajectories" for three generic slow 
roll models, \<j) A (dotted), m 2 (j) 2 (solid), and a representative natural inflation model 
(dashed), along with the position at different values of N — 30 (triangle), N = 40 
(square), N — 50 (star), N = 60 (circle). 
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Figure 4. The joint 68% (dark) and 95% (light) confidence levels obtained on the 
first three HSR parameters {e, 77, ^} at the fiducial scale fc* = 0.02 Mpc^ 1 from the 
WMAP 5 year data, assuming all higher order slow roll parameters are zero. The 
grey (dotted) contours correspond to a "minimal" e-fold prior, N(k+) > 15. The blue 
(solid) contours are obtained by assuming instantaneous reheating. 
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Figure 5. The joint 68% (inner) and 95% (outer) bounds on the first three HSR 
parameters {e,»7,C} a t the fiducial scale fc* = 0.02 Mpc" 1 . The blue constraints 
are derived from WMAP 5 year data alone, and the red constraints from the 
WMAP5+SNLS data combination, with T leh > 10 TeV. 
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Figure 6. The joint 68% (inner) and 95% (outer) bounds on the power law spectral 
parameters at the fiducial scale fc* = 0.02 Mpc -1 , obtained by transforming the 
constraints shown in Fig.[5]into this parameter space. The blue constraints are derived 
from WMAP 5 year data alone, and the red constraints from the WMAP5+SNLS data 
combination. The dotted contours come from Slow Roll Reconstruction applying an 
e-fold prior assuming T rc h > 10 TeV. For comparison, the dashed contours show an 
analysis using the empirical power-law prescription in terms of n s , r, and dn s /dh\k 
at a pivot scale of 0.02 Mpc" 1 , using WMAP5 data. 
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4. Primordial Black Holes and the Risk of Eternal Inflation 

In contrast to the constraints found with the WMAP3 dataset [El [7J, the WMAP5 
results [25J show less evidence for a strong running of the scalar spectrum. While 
second order terms in slow roll contribute to the running, these cannot be resolved with 
current astrophysical data [BJ. Therefore the running is dominated by the third slow 
roll parameter £, and in this limit, 

dn* , , 

a ^^= 3 - 2 «- < 13 » 

Large, positive values of £ lead to a very short period of inflation, and are thus excluded 
by an e-fold prior [7J and, as with WMAP3, the upper limit on £ derived from WMAP5 
is driven solely by the e-fold prior. The lower limit on £ from WMAP5 is actually more 
negative than the analogous result from WMAP3 - so the maximal permitted running 
is thus more positive. The difference is primarily due to a better treatment of the beam 
and point sources in the latest release, rather than the extra two years of data. This 
situation illustrates the importance of systematic effects, and reminds us of the need for 
caution when interpreting cosmological constraints. 

Interestingly, significantly negative values of £ can also be excluded on physical 
grounds, again under the assumption that the first three HSR parameters give a 
description of the inflaton dynamics till the end of inflation. With a substantial negative 
£, we find a class of solutions with e — > 0, and the field rolls towards a minimum with a 
substantial vacuum energy. Given that e — > while H remains finite, the amplitude of 
the perturbation spectrum diverges. As Pr approaches unity (at small scales) primordial 
black holes can form in the post-inflationary universe [2H1 EHl EH EH EH [321 E31 EU [35] . 
Their mass is limited by the energy contained within the horizon volume, which increases 
as the universe expands. A black hole formed at a temperature greater than ~ 10 8 GeV 
decays before the present epoch. Consequently, if stable black holes form after inflation 
with a high (near-GUT) reheat temperature, they are generated by modes which reenter 
the horizon with a substantial amplitude after the temperature has fallen below this 
critical value. In particular, [32J, [33J, [35] consider primordial black hole formation after 
slow roll inflation, and the impact of a running scalar index driven by a large negative £ 
term. Interestingly, [33] uses a version of Monte Carlo reconstruction [101 [11] to look for 
potentials whose perturbations have spectral indices consistent with CMB data while 
also producing large numbers of small black holes. The only potentials they find which 
fit both criteria have sharp "features", and thus do not fit naturally within the slow 
roll hierarchy. The authors of [35] note that this conclusion is at odds with|32] . and 
attribute the discrepancy to "the use of flow equations and the hierarchy"!^ by [33J. 

% By "the hierarchy" , which in standard usage means the infinite hierarchy of differential equations 
describing the evolution of the HSR parameters, the authors of |35j instead mean progressively tighter 
restrictions on the initial conditions for the HSR parameters as given by their equation (55). However, 
[33j only apply these restrictions to the priors on their initial conditions. Individual draws from their 
priors are manifestly able to generate dynamics which violate these restrictions. Thus the objections 
of [35] do not apply. 
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However, as we saw in Figures and El our chains - which also use the flow equations - 
have no difficulty finding regions of parameter space which lead to the overproduction 
of primordial black holes. 

Instead, the discrepancy is attributable to the way that "Monte Carlo 
reconstruction" [TQl E] handles models in which inflation can continue indefinitely. 
Namely, one assumes that if there is a late time attractor, the inflationary trajectory 
reaches this attractor, whereupon inflation ends via a hybrid transition. As noted in 
[TT] . this is a simplifying assumption. It ensures that the spectrum has a blue spectral 
index at CMB scales, since (in slow roll terminology) V is vanishingly small and V" 
is negative, and these models are excluded by the data. Thus, one cannot find sets of 
slow roll parameters which overproduce primordial black holes while also matching the 
CMB power spectrum with this ansatz. Conversely, Slow Roll Reconstruction works 
with the values of the slow roll parameters at CMB scales, and effectively marginalizes 
over the mechanism by which inflation ends. We thus recover models which have not 
quite reached the attractor solution e — > at CMB scales, so their spectral index is 
red with a large positive running. By the end of inflation, the spectra generated by 
these models are apparently large enough to source the production of primordial black 
holes, and this class of model is missed by the algorithm used by [33] . In other words, 
they make an arbitrary and non-unique choice for the piece of the trajectory which 
corresponds to the CMB scales which a priori rejects the models we have found which 
both match constraints at CMB scales and lead to primordial black hole overproduction 
by the end of inflation. The running mass model is an explicit inflationary scenario with 
these dynamics [55l 156] , invoking softly-broken global supersymmetry during inflation. 
This model has been compared to observations in a number of works [53, [58] [32] [59], 60J . 

In what follows, we examine the lower limit on £ that results from imposing an 
additional prior based on primordial black hole overproduction. The limit on the 
primordial amplitude at the end of inflation is somewhat uncertain - a variety of 
arguments [61] [62] [30] [63] EU [35] suggest that P-r < 10~ 3 -10~ 2 . Further, as we have 
seen, in the absence of a concrete mechanism to end inflation, there is a large uncertainty 
in the identification of where inflation ends in these models. As discussed previously 
(for example, by Leach and collaborators [32J ) , this uncertainty has a significant effect 
on the constraint as well. To take into account these uncertainties, we post-process our 
chains as follows. We identify the end of inflation N end using equation (TTTj) assuming a 
minimal T re h = 10 MeV and a maximal requirement of instantaneous reheating. Further, 
at N end , we apply two amplitude cuts: Pn < 1CT 3 and P^ < 1CT 2 . As pointed out in 
[32] , the Stewart-Lyth formulae are accurate enough to compute the power spectrum at 
the end of inflation in these models, provided that the background solution is computed 
exactly, as we have done. We show our results in Figure [7] We see that the lower 
limit on £ is strongly dependent on the precise form of the prior. For all but the weakest 
combinations of priors, a significant part of the observationally-allowed parameter space 
is excluded by primordial blackhole overproduction. 
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Figure 7. The joint 68% (inner) and 95% (outer) bounds at the fiducial scale k+ = 
0.02 Mpc -1 . The red constraints come from the WMAP5+SNLS data combination 
using Slow Roll Reconstruction, applying an e-fold prior assuming T rc h > 10 TeV. 
This e-fold prior eliminates models with large positive £ and has no effect on models 
with large negative £, which inflate forever in the absence of a hybrid transition to end 
inflation. The dark contours come from post-processing the same chains to exclude 
models which lead to primordial black hole overproduction (see text for discussion): 
solid black contours assume that T re h = 10 MeV, and the dot-dash contours assume 
instantaneous reheating. The left and right hand columns are based on requiring 
that the power spectrum amplitude at the end of inflation satisfy Pr. < 10~ 3 and 
P-ji < 10 -2 respectively. We thus see that the black hole constraint depends strongly 
on the duration of inflation and the precise value of the amplitude at which excessive 
black hole production begins. 



Even if we ignore the previous constraint (though it should be clear that the 
following case is subset of the former), eternal inflation occurs in models where the 
quantum fluctuations (which are proportional to H) dominate the classical rolling, and 
this can be shown to be equivalent to the density perturbations being of order unity 
(36j [35] . In this case we could not trust the usual classical inflaton equations of motion, 
which are embodied in the slow-roll hierarchy, and inflation cannot end coherently. 
This again rules out large, negative values of £ in three parameter models. The eternal 
inflation prior is implemented as follows. First, we note that in the models where there 
is a danger of eternal inflation, e — > while H(<f>) asymptotes to a constant, i^asymp- As 
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Figure 8. The joint 68% (inner) and 95% (outer) bounds at the fiducial scale k+ = 
0.02 Mpc -1 . The red constraints come from the WMAP5+SNLS data combination 
using Slow Roll Reconstruction, applying an e-fold prior assuming T rc h > 10 TeV. The 
solid black contours come from post-processing the equivalent WMAP5-only chains to 
impose an additional prior based on avoiding eternal inflation (see text for discussion) . 
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we approach this limit, the primordial power spectrum H 2 /e diverges. We assume that 
H cnd = H reil and H cnd = -f^asymp, where the latter quantity can be computed analytically 
for the {e,77,£} parameterization by solving a quadratic equation to find the zero of e. 
Then we substitute this H end into equation (TTTT) to obtain N end . The eternal inflation 
prior is then implemented by requiring that -P"R.(iV en d) < 1. Figure M shows the effect 
of applying this prior. We see that the WMAP5+SNLS constraints with no eternal 
inflation prior are roughly the same as the WMAP5 data combination with this prior 
- i.e., WMAP5+SNLS is able to rule out the eternally- inflating part of the parameter 
space which is allowed by WMAP5 alone. 

5. Discussion 

We present constraints on the slow roll parameters, following the WMAP5 data release 
[221 E3l El [25], in combination with ACBAR (ACBAR) [26], and SNLS [27]. This 
analysis updates our treatment of the WMAP3 [51 [6] release. In addition, we have 
presented a more careful implementation of the "e-fold prior" , which we use to eliminate 
models that do not produce sufficient inflation. Further, we have added a prior on the 
inflationary parameter space that excludes models which over-produce primordial black 
holes, or would permit the onset of eternal inflation after CMB scales have left the 
horizon. 

We can easily post-process our chains to present our results in terms of the usual 
{n s , r, dn s / d In k} formalism. With two slow roll parameters, the permitted ranges of n s 
and r closely match those obtained by including these parameters directly in the chains. 
Conversely, with three slow roll parameters, Slow Roll Reconstruction leads to very tight 
limits on dn s /d\nk - especially when we include a primordial black hole constraint. It 
is important to realize that these limits are being driven by our priors, rather than the 
data. The ability of Slow Roll Reconstruction to include these priors - especially when 
they relate to the duration of inflation - is a key advantage of this methodology. On 
the other hand, despite the great strides that have been made in recent years, current 
cosmological data does not significantly constrain £, or V", given that the bounds on 
this parameter are largely functions of the priors. We note that the precise form of the 
primordial black hole bound is currently somewhat ambiguous, and a thorough analysis 
of this constraint should be undertaken in the future. However, even with current data, 
the combination of WMAP5 and SNLS essentially eliminates the region of parameter 
space that would support eternal inflation. Likewise, with weak assumptions about 
reheating, the region of parameter space which leads to the overproduction of black 
holes is roughly similar to that excluded by WMAP5+SNLS. However, if we assume the 
reheating scale is close to the GUT scale, the prior rules out regions of parameter space 
currently allowed by the data. 

One can view the overlap between the region of the {n s , dn s /d\nk} plane permitted 
by a conventional analysis and Slow Roll Reconstruction shown in the lower panel 
of Figure | as a consistency check - if these two regions did not overlap, we could 
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exclude all 3-parameter models of single field inflation. Both eternal inflation and the 
overproduction of primordial black holes can be avoided by including a fourth slow roll 
parameter, which makes £ scale dependent. We do not explore this option here: since £ 
is only weakly constrained by the data, the bounds on a fourth slow roll parameter ( 3 A) 
would be almost entirely set by the priors. Further, given that these parameters are 
not strongly constrained by the data, the central values we compute will be functions 
of our prior distributions. Since e effectively measures the inflationary scale (once the 
amplitude of the scalar perturbation spectrum is known), a flat prior on e biases us 
toward high inflationary scales, whereas a log prior would tip us towards inflationary 
scales that are that are below the GUT scale - with the lower bound coming from our 
assumptions about reheating, as explained in Section 2. Conversely, different slow roll 
parameterizations can suggest a "detection" of tensor models within current data [65J. 
This ambiguity is a reminder of the importance of priors in Bayesian statistics [66J and 
is characteristic of a parameter which is not (currently) well-constrained by data. 

Needless to say, we expect this situation to improve in the future, but given that 
the permitted region of the {n s , dn s /d In k} plane is currently much larger than that 
consistent with the slow roll prior, we need more than incremental improvements in the 
global cosmo logical dataset. It is clear that Planck will lead to a substantial tightening 
of constraints on the running [671 [68] . Further, if dedicated polarization experiments 
detect primordial tensor modes this would put a lower bound on e, removing most of the 
ambiguity in the inflationary scale, and thus the e-fold and reheating priors built into 
Slow Roll Reconstruction. Moreover, we can expect similar improvements from large 
scale structure measurements, supernovae, or measurements of the primordial spectrum 
from high redshift 21 cm experiments. In particular, if e is not vanishingly small, second 
order terms in slow roll naturally lead to a small running of the scalar spectral index. 
This running cannot be detected in CMB data alone [Sj, but may well be detectable 
in combinations of several datasets, providing further constraints on the inflationary 
parameter space. 
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Appendix A. Parameter Constraints from Monte-Carlo Markov Chain 
Analysis 
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Table Al. Constraints on HSR parameters defined at k = 0.02 Mpc" 1 for the 
{e, 77, log[10 10 ^4 s ]} parameterization. Constraints are at the 68% confidence level unless 
otherwise noted. 
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< 0.015 (95% CL) 


V 


ni7+o-022 

U - Ui ' -0.022 


012+ - 019 


006+ - 017 

U.UUO_ 016 




-0 nn7+ - 006 

U.UU 1 _ .006 


-o no c i +0 - 005 

U.UUO_o.oo5 


-0 00fi +a005 

U.UUU_ 005 


ln[10 10 A s ] 


o nc:+0.05 
°- uo -0.05 


o Q7+0.04 
°- u ' -0.04 


o Q7+0.04 
°- u ' -0.04 



Table A2. Constraints on HSR parameters denned at k = 0.02 Mpc^ 1 for the 
{e, f], £, log[10 10 A s ]} parameterization. Constraints are at the 68% confidence level 
unless otherwise noted. 



